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Fractal geometry of critical curves appearing in 2D critical systems is characterized by their 
harmonic measure. For systems described by conformal field theories with central charge c ^ f , 
scaling exponents of the harmonic measure have been computed by B. Duplantier [Phys. Rev. Lett. 
84, 1363 (2000)] by relating the problem to boundary two-dimensional gravity. We present a simple 
argument connecting the harmonic measure of critical curves to operators obtained by fusion of 
primary fields, and compute characteristics of the fractal geometry by means of regular methods of 
conformal field theory. The method is not limited to theories with c ^5 1. 
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Introduction. Two-dimensional (2D) statistical mod- 
els typically exhibit stochastic curves, such as external 
perimeters of critical clusters in the Potts model or fluc- 
tuating loops in models of Refs. [J-Q- In the critical 
regime these curves are fractal. Different phases of loop 
models exhibit two distinct types of critical curves: di- 
lute and dense. Dilute curves are simple while dense 
curves have infinitely many double points at which they 
touch but do not intersect. There exists a duality re- 
lation between the two types, in particular, the external 
perimeter of a dense curve is a dilute curve [J] . Despite of 
a long history of studying critical phenomena, inquiries 
into the stochastic geometry of the critical curves are rel- 
atively recent. Traditional conformal field theory (CFT) 
concentrates on correlations of local fields and essen- 
tially leaves their relation to geometrical objects, such 
as critical curves, unclear. B. Duplantier, in a seminal 
paper found the harmonic multifractal spectrum of 
critical curves arising in critical systems with the cen- 
tral charge c ^ 1 (see Eqs. dH [8])). His method used 
an intriguing connection between CFT and 2D boundary 
quantum gravity. More recently, the stochastic Loewner 
evolution approach was successfully used for the same 
purpose 0. Both methods are interesting and powerful, 
but somewhat foreign to traditional CFT approach and 
not obviously susceptible to generalizations. 

In this letter we show that geometrical properties of 
dilute critical curves are naturally linked to correlation 
functions of primary fields and can be studied within the 
regular framework of CFT, not limited to theories with 
c ^ 1. Our arguments may also shed some light on the 
boundary quantum gravity itself. 

Harmonic measure of critical curves. A basic char- 
acteristic of simple curves is the harmonic measure ^. 
It has a simple electrostatic interpretation. Consider a 
closed curve 7 made of a conducting material and carry- 
ing a total electric charge of one. Harmonic measure of 
any part of the curve is the charge of this part. In what 
follows we will pick a point of interest zq on the curve 
and consider a disc of a small radius r centered at zq as 
is in the Figure. It surrounds a small part of the curve, 
and we define fi{r) to be harmonic measure of this part. 

If a curve is a closed loop lying entirely in the bulk (i.e. 



it does not touch the system boundaries), all its points 
are statistically identical. We can also consider the case 
when a curve emerges from a boundary, and the case 
when a curve has one or both dangling ends in the bulk. 
The latter is achieved by inserting into the system a small 
isolated boundary (almost a puncture) so that the curve 
ends on this boundary. Statistics of harmonic measure 
is different in the following three cases Q depicted in 
the Figure: (i) bulk, if zq lies in the bulk and is not 
an endpoint; (ii) boundary, if Zq is a point where 7 is 
connected with the system boundary; (iii) extremity, if zq 
is a dangling endpoint of 7 in the bulk. More generally 
one may consider n curves emanating from the boundary, 
or n curves meeting at a point in the bulk. The bulk and 
extremity cases correspond to n = 2 and n = 1. 
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FIG. 1: Harmonic measure of a critical curve evaluated in 
the bulk, on the system boundary, and at an extremity. 



Let us define a conformal map w{z) of the exterior of 7 
to some standard domain. For a closed loop in the bulk it 
can be the exterior of a unit disc, while for a curve with 
both ends on the system boundary or a "sprout" , the 
upper half plane is more convenient. We normalize the 
map so that the point of interest zq is mapped onto itself, 
choose it to be the origin zq = 0, w{0) = 0, and demand 
that w'{oo) = 1. The moments of harmonic measure 
of dilute critical curves scale at small r with non-trivial 
exponents. The relations fi{r) ^ \w{r)\ ^ r\w'{r)\ follow 
from the definition of /x(r). Therefore, the scaling of the 
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moments of fi{r) is determined by that of [^'(r)!: 

^|u;'(r)|V- r^(''\ r ^ 0. (1) 

The notation -<. . stands for statistical average over the 
ensemble of critical cm'ves. The notation (...) is reserved 
for correlators of CFT. 

Harmonic measure and fluctuating geometry. The 
scaling of the harmonic measure may be computed by 
considering CFT correlation functions. It is the easiest 
to start with a point where a curve 7 connects with the 
system boundary. We assume that the system occupies 
the interior of an annulus. The curve starts from the 
large circle at zq = 0, ending on the large circle at a 
point L. 

The partition function Z{0, L) restricted to configura- 
tions that contain a curve 7 connecting the points and 
L, is given by the correlator of two appropriate boundary 
operators B changing boundary conditions Q . These op- 
erators create curves emanating from the boundary [lO|: 



Z{Q,L)/Z={B{0)B{L)), 



(2) 



where Z is the unrestricted partition function. The sub- 
script of (...) refers to the domain of definition, in this 
case the upper half plane. 

This correlator, as indeed any correlator containing 
two boundary operators B, can be computed in two steps: 
In the first step we pick a particular realization of the 
curve 7. Within each realization, the curve 7 is the 
boundary separating two independent systems — the inte- 
rior and the exterior of 7 — with partition functions Z"'* 
and respectively. These are stochastic objects that 
depend on the fluctuating geometry of 7. At the second 
step we sum over the realizations of 7. We thus obtain 
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(3) 



We further insert an additional boundary primary op- 
erator Oh{r) of conformal weight h sufficiently close to 
0, and a second copy Oh{oo). Both act as a source of 
harmonic measure. We thus consider the correlator 



(O^(r)S(0)6(L)O^(oo))i 



(4) 



Since we are only interested in the j'-depcndcnce of the 
correlator, we can fuse together the distant primary 
fields: B{L) x 0^(00) = \E'(oo). We therefore consider 
the r-dependence of a 3-point function 



(O,,(r)6(0)*(oo))H, 



(5) 



and show that it yields the statistics of the harmonic 
measure. 

Decomposing the upper half plane into the exterior 
and the interior of 7 as before, we can rewrite (|4]) as the 
average over the fluctuating geometry of 7: 



Here the domain of the definition of the correlator of 
primary fields is the exterior of 7. This correlator is sta- 
tistically independent from the other two factors in the 
numerator of ([6]) in the limit r <C \L\. We are left with 
the correlation function {Oh{r)Oh{oo))^^ of two primary 
fields of boundary CFT, further averaged over all config- 
uration of the boundary 7. It is equal to the 3-point 
correlation function ([5|). 

Now we apply the conformal transformation w{z) 
which maps the exterior of 7 onto the upper half 
plane. Being a primary operator of conformal weight 
h, Oh{r) transforms as Oh ^ Oh{w{r))\w' {r)\'^ , while 
Oh{oo) does not change because of the normalization 
of w{z). The transformation relates the correlation 
function in the exterior of 7 to a correlation func- 
tion in the upper half plane: {Oh{'w{r)) Oh {00))°^^ = 
\w' {r)\'^ (Oh{w{r))Oh{oo))M- The latter does not depend 
on r, but rather on the size of the entire system. 

Summing up, we obtain a scaling relation between the 
moments of the harmonic measure near the boundary 
and correlation functions of primary operators |ll| : 



(O,,(r)6(0)*((X3))H ^ ^|u.'(r)| V, r « (7) 

This scaling relation is the main result of the paper. It al- 
lows us to reproduce the scaling of the harmonic measure 
upon identification of the operators B. 

The spectrum of the harmonic measure for CFTs with 
c ^ 1. The scaling exponents of -< |zi;'(r)|'' >- in the 
bulk, near the boundary and near an extremity are de- 
noted AbuikCi), ^{h) and Aoxtr(^), respectively. If we 
parametrize the central charge of the critical system by 
c = 1 — 6(\/4/k — where < k < 4 for dilute 

and K > 4 for dense critical curves, then the exponents 
obtained in [3| read 

Acxtr = -V2 + Ak/8, Abulk = Acxtr + a/2 (8) 

The boundary exponent happens to be identical to the 
dressed gravitational dimension. This is the dimension 
of a primary field on a random surface, which dimension 
on a flat space is h. They are connected by equation Q 
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(9) 



<{Oh{r)Oh{oo)YfZ'j^''Z';^V IZ 



(6) 



In order to derive these results along the lines presented 
above, we must identify the operators creating critical 
curves. In what follows this result and its generalizations 
to the case of n curves meeting at the same point, will 
appear in a compact form in Eqs. (|15l lISp in terms of 
the charges of relevant operators. 

Levels of a Gaussian field. The relation between criti- 
cal curves and boundary CFT is most transparent in the 
case when they are dense (k > 4) upon representation 
as level lines of a Gaussian field Below we 

discuss this representation and find the operators rep- 
resenting dense critical curves. The duality then implies 



3 



that the corresponding operators for the dilute curves are 
obtained by simple continuation to k ^ 4. 

A c ^ 1 CFT in an annulus D is described by a 
real compactified Bose field ip{z,z). In the normaliza- 
tion where the radius of compactification is equal to 1, 
so that ip ~ tp + 2n, the classical action reads [3 : 

Here K is the geodesic curvature of each component of 
the system boundary and the stiffness g and the "charge" 
Co are related to k as f/ = 4/k, cq = 1 — 4/k. The last 
term in the action (jlOp represents the marginal part of a 
27r-periodic locking potential. 

We impose the boundary condition 

dnP\dD = 0, (11) 

where the derivative is normal to the boundary. The real 
field Lp is the sum of a holomorphic and an antiholomor- 
phic parts and a real zero-mode (j)Q : 

p{z,z) = (j){z)+J{z) + <i)a. (12) 

The fields 0(2;) and 4>{z) arc glued together on the bound- 
ary by the condition (fTTj) which reads di(j){z) = di(j){z), 
where di is tangential derivative. Both parts can be 
considered as one holomorphic field on a torus (Schot- 
tky double) obtained by gluing the annulus to its re- 
flected copy along the boundaries. The "dual" Bose field 
'p{z,z) = {A/iK){4i{z) — ^(z)) is related to p> through 
Cauchy-Riemann conditions. 

The gradient J = Vp){z, z) is therefore the current. 
It is conserved. On the boundary [zj = 1 the boundary 
condition (|lip then means that no vector current fiows 
through the system boundary Jn\dD = 0. 

Fluctuating loops. The critical curves have a sim- 
ple interpretation in terms of the Bose field. They are 
the level lines Re(^(z,z) = kn of the height function 
Re(/3(z,z). The level lines are non-intersecting plane 
loops which are identified with boundaries of critical clus- 
ters in the Potts model, or lines in the 0{n) model 
01 . In this formulation, statistical models represented by 
CFTs with c ^ 1 can be seen as a gas of fluctuating loops. 
The boundary condition ensures that the loops cannot 
cross the system boundaries and, furthermore, that the 
system boundaries themselves are loops. 

In the Hamiltonian formalism of radial quantization, 
the non-contractible loops C in the annulus represent 
coherent states propagating along the cylinder C = 
logz. Coherent states are defined by the condition 
tp(z , z)\C , 4>q) = 0o|C,0o)i z eC. Contractible loops rep- 
resent virtual states. The partition function with the 
action pO|) can be seen as the overlap of the boundary 
coherent states. 

A normalization of </3, customary in the CFT litera- 
ture, fixes 5 = 5 and introduces the notation oij^ — 



•\/ 4/k, Q!_ = — ■\/k/4, 2q;o = a+ -l- a- 0]. The radius 
of compactification is then TZ = ^JS/n = \^a+. Below 
we proceed in the physical normalization, where TZ = I. 

Primary operators. In terms of the Bose field the pri- 
mary operators read 

0('='")(z, z) = e*«v(^^^)e™'^(^'^), (13) 

where e and m are "electric" and "magnetic" charges. 

The holomorphic weight of this operator is h{a) = 
o? — 2aaoj where a = — i(eQ;_ -I- TOa_|_) is the holo- 
morphic charge. The antiholomorphic charge of is 
a = — i(eQ;_ — toq;+). It is also customary to label the 
primary operators and their holomorphic charges by two 
numbers r, s using the Kac table Q!r,s = ^(1 — + 
i(l — s)a_. Here r = l-|-m, s = l + e. 

The magnetic operator with e = 0, jti = 1 introduces 
a defect line on which the value of the field p changes 
by the compactification length 27r. The electric operator 
with e = 1, m = picks up a phase difference 27r while 
going around the magnetic operator. A general C^"^'™) is 
the composition of the two. 

Operators representing critical curves Q, 0, [7^ /. The 
magnetic operator with charge m, applied at the bound- 
ary, changes the values of tp by Trm (since one can go 
from one side of the operator to the other by half of a 
fuU turn): O(0''")(0)(/?(a;) = ((^(x) -f m7r6'(a;))C'(°'") (0), 
where 9{x) is the step function and x is the coordinate 
along the boundary. Thus, it is a boundary condition 
changing operator. In particular, 0^^'^\Q) changes the 
boundary condition by half the compactification length, 
just as a cluster wall does, so it produces a critical curve 
emanating from the system boundary at x = in the 
direction normal to the boundary. This is, therefore, the 
previously mentioned operator B. It happens to be de- 
generate on level 2 and is placed in the Kac table as[3| 
'I'2.i- We can insert two such operators on the external 
boundary of the annulus and the other two on the cen- 
tral puncture in order to connect the boundaries with 
two curves. These curves then divide the annulus into 
two domains, each having the topology of a disc. 

Similarly, the boundary operator O^*^-") = "^^^^i i is 
degenerate on level n -I- 1 and produces n curves. Its 
holomorphic charge is a„+i,i = — ^q;+. 

Pinning operator. This operator creates n curves em- 
anating from a point in the bulk, that is, from a punc- 
ture. The puncture, as any internal system boundary, 
carries the electric charge — ep. The combination of this 
charge and the magnetic charge associated with the cre- 
ation of n curves identifies the bulk pinning operator 
as 0(~'=O'"/2) with the conformal charge ai4.„/2,4/K = 
oi^fi = —ja^ + aQ. In the Kac classification it is the 
fleld ^-^,0. 

Screening operators. The bulk magnetic operator 
with the double magnetic charge O^^'^^ cuts a loop and 
changes the direction of each part by ±7r giving rise to a 
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new loop representing a virtual state. This operator must 
be marginal, which is the origin of the relation between 
the stiffness and the screening charge eo+g = 1. Another 
marginal (i.e. screening) operator 0(2,0) represents the 
marginal part of the locking potential in (jlOp . 

Boundary exponents. We argued that the curve- 
creating operator ;B in ([7]) is ^'2,1. Therefore, ^ reads 

{Ohir)^2,m'^i^))a ~ ^k'WI V, (14) 

where the holomorphic charge of Oh is Uh = ao — 
y^a^+h and that of 'J (00) is 2ao — ah — 012,1 (the choice 
of dual charges is made unique by A(0) = 0.) 

On the other hand, the scaling behavior of this cor- 
relator is easily computed by regular means of Coulomb 
gas technique It scales as 7-20^02,1^ ^j-^g comparison 
yields the result ^ written in a suggestive form: 

A(/i) = 2a/ia2,i, ah = a^- \J al + h (15) 

where string susceptibility 7string = 1 ^ 4q;2 1 . 

An immediate generalization of this formula A^") = 
2ahOin+i,i — «A can be obtained by replacing '^2,1 by 
in (|14p . It describes the harmonic measure of n 
curves reaching the system boundary at one point. 

Bulk exponents. Similar arguments yield the value of 
the bulk exponents. Let points and L lie in the bulk 
and consider a bulk correlator 

(0,,.(r)*i,o(0)*i,o(L)0,,.(oo)). (16) 

The fields ^'i_o ensure the existence of a closed curve 
7 connecting and L. In the limit r <C \L\ and be- 
cause we are only interested in the r-behavior, we can 
fuse '^ifl{L) X Oft' (00) = 4'(oo) where the charge of ^ is 
2q!o — Oih' — Oiifi, and consider instead a 3-point function 
(Oft'(r)«'i,o(0)'«'(oo)). 

As in the boundary case, we argue that this correla- 
tor equals <{Oh' {r)Oh' {oo))'^'^>- up to a normalization. 
We transform it into a correlator in a disc exterior via a 
map w{z) which takes 7 to a unit circle centered at — i 
such that Oh'ir) \w' {r)\'^^' Oh'{w{r)). The difference 
with the boundary case is that Oh' is now a bulk field in 
the presence of a circular boundary. We therefore fuse 
Oh'{w{r)) with its image inside the disc (approximately 
at w*{r)) and take the leading non-trivial fusion product 
{'w{r) — ii;*(r))''^2'i 0^(0). Here Oh is a primary field of 
weight h. We determine h through the neutrality condi- 
tion 2ah' + (2ao - ah) = 2ao, yielding ah = 2ah'. 

Since r is small, 'w{r) — w*{r) ^ r\w'{r)\. Summing up, 
we give the scaling behavior of the original correlator: 

(Ov(r)*i,o(0)vl/(oo)) ~ r''-2h' ^|^'(^)|V . (17) 

On the other hand, scaling behavior of this correlator at 
small r is easily found by regular CFT means: it scales 
as 7-4Qh/ai,o^ -y^g ^Yius obtain 

Abuik(/i) = 2/i'-/i-H4aft'ai,o, (18) 



which coincides with the result given in Eq. ([8]). 

A replacement XPi^o —^ ^^,0 in (HH) produces n curves 
emanating from a given bulk point. Their scaling expo- 
nents are obtained by replacing ai^o "^f ,0 in 

^iti^W = -h/2 + aft(2a„/2,o - ao)- (19) 

If n is even it gives a bulk exponent for the case when 
n/2 curves are passing through one point. The case n = 1 
yields the extremity exponent Acxtr(^) of a single dan- 
gling end in the bulk 

Discussion The method of computing the statistics 
of the harmonic measure, discussed here, is amenable 
to generalizations. For instance, one can compute mul- 
tipoint correlation functions. Another generalization is 
for curves generated by conformal field theories with 
c > 1 [3, 13]. In particular, it has been shown in 
[ist that SU(2)fe - Wess-Zumino-Witten model generates 
critical curves identical to those generated by its coset 
su(2)fc©su(2)i/su(2)fc+i — a unitary minimal model with 
K = 4|i|. Its substitution into Eqs. ([9l[8]) gives the scal- 
ing of the harmonic measure for this case. 
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